Spin Hall effect in the kagome lattice with Rashba spin-orbit interaction 
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We study the spin Hall effect in the kagome lattice with Rashba spin-orbit coupling. The conserved 
spin Hall conductance a^y (see text) and its two components, i.e., the conventional term crj^ and 
the spin-torque-dipole term crj^, are numerically calculated, which show a series of plateaus as a 
function of the electron Fermi energy ep. A consistent two-band analysis, as well as a Berry-phase 
interpretation, is also given. We show that these plateaus are a consequence of the various Fermi- 
surface topologies when tuning ep. In particular, we predict that compared to the case with the 
Fermi surface encircling the F point in the Brillouin zone, the amplitude of the spin Hall conductance 
with the Fermi surface encircling the K points is twice enhanced, which makes it highly meaningful 
in the future to systematically carry out studies of the K-valley spintronics. 
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Spintronics, which combines the basic quantum me- 
chanics of coherent spin dynamics and technological ap- 
phcations in information processing and storage devices, 
has been become a very active and promising fieldiii^i^. 
The key is how to control and manipulate the spin de- 
grees of freedom. One of the tools is using the spin- 
orbit (SO) couphngs, which describe the interactions be- 
tween the electron's orbital and spin degrees and pro- 
vide an ability to manipulate the spin state via chang- 
ing some external factors, such as an external electric 
field. It has been argued that the SO interaction leads 
to an intrinsic spin Hall effect (SHE}^i^, in which a 
spin current flows perpendicular to an applied electric 
field. The initial theoretical^i^iii^ia>iaiiiii^ii^»ii^ and 
experimenta l ""^^i^^'^^i^^ studies of SHE were mainly fo- 
cused on the p OT n doped semiconductors (such as 
GaAs). Then, Murakami et. al^^ first identified a class of 
cubic materials that are usual insulators, but nonetheless 
exhibit a finite spin Hall conductance (SHC). In those 
proposed "spin Hall insulators" (SHIs) the SHC is not 
quantized and depends on the system parameters. Later 
and even more fundamentally, it has been evolving into 
one important theme in condensed matter physics that 
the SHC can be quantized in time-reversal invariant sys- 
tems and thus can be used as an order parameter to char- 
acterize the emergence of new topological insulating state 
of mattc r^^-^^-^'^-^'^-^^-^^-^''-^^'^^-^"''^^-'^^''^'^-^^-'^^ 

It is clear now that besides the external SO coupling 
(e.g., the Rashba SO couphng), the lattice structure itself 
also has crucial impact on SHE through the related band 
structure. Different lattice structure may produce new 
features in the spin transport, which provides versatile 
choices of materials to study spin Hall transport. Mo- 
tivated by this observation, in this paper we study the 
intrinsic SHE of the noninteracting electrons in a two- 
dimensional (2D) kagome lattice with Rashba SO cou- 
pling. Since our attention is solely on the SHE character 
brought about by the interplay between the kagome lat- 
tice structure and the Rashba SO coupling, thus unlike 



most of previous works, the kagome lattice considered in 
this paper is nonmagnetic. The nonmagnetic kagome lat- 
tice structure has been either fabricated by modern pat- 
terning techniques^i^ or observed in reconstructed semi- 
conductor surfaces^. In the former case, remarkably, the 
electron filling factor (namely, the Fermi energy) can be 
readily controlled by applying a gate voltage^^. Our lat- 
tice model is free from the constraint imposed on the 
k-p approximation used in the extensively studied GaAs 
two-dimensional electron gas (2DEG), in which the k-p 
Hamiltonian is only valid around the F point in the Bril- 
louin zone (BZ). In contrast, our lattice model allows for 
any electron filling, which result in various Fermi-surface 
topologies, which in turn, as will be shown below, pro- 
duces profound effects on the spin Hall transport. 

To calculate the SHC and build a correspondence be- 
tween spin current and spin accumulation in the present 
SO coupled system, in which the electron spin (sz here, to 
be specific) is not conserved, we use a "conserved" spin 
current J/r^, which is a sum of the conventional spin 
current Js=i{v, s^} and a spin torque dipole PT-=rs2. 
This spin current satisfies both spin continuity equation 
i9tSz+V- Js=0 (within spin relaxation time) and Onsager 
relation^^. If the spin itself is conserved (as in quantum 
SHIs) , J^s is reduced to . In general, the spin transport 
coefficient a"^^ under new definition is composed of two 
parts, i.e., the conventional part cr^^ and the spin torque 
dipole correction cr^l^. A general Kubo formula^ for the 
spin transport coefficients is employed in this paper to 
calculate the SHC. 

Let us consider the tight-binding model for indepen- 
dent electrons on the 2D kagome lattice (Fig. 1). The 
spin-independent part of the Hamiltonian is given by 

i^o=io5I(cLcja+H.c.), (1) 

where tij=to is the hopping amplitude between the near- 
est neighbor link c]^ {cia) is the creation (annihi- 
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FIG. 1: (Color online) Schematic picture of the 2D kagome 
lattice. The dashed lines represent the Wigner-Seitz unit cell, 
which contains three independent sites (A, B, C). 



lation) operator of an electron with spin a (up or down) 
on lattice site i. For simplicity, we choose to=l as the 
energy unit and the distance a between the nearest sites 
as the length unit throughout this paper. 

The Hamiltonian ([T|) can be diagonalized in the mo- 
mentum space as 



^0 = V7/^l[(Ho(k)0l2x2)V'k, 



(2) 



where the 2x2 
spin degeneracy 



unit matrix I2x2 denotes the 
in the Hamiltonian Hq, and 



V'k=(c^kT,CBkT,cckT,CAki,cski,ccki) is the six- 
component electron field operator, which includes the 
three lattice sites s {=A, B, C) in the Wigner-Seitz unit 
cell shown in Fig. 1. Each component of V'k is the 
Fourier transform of Qq,, i.e.. 



■080 (k) 



mn 



(3) 



where we have changed notation i-^{mns) by using (mn) 
to label the kagome unit cells, -ffo(k) is a 3x3 spinless 
matrix given by 

(0 2cos(k-ai) 2cos(k-a3) 

2cos(k-ai) 2cos(k-a2) 

2cos(k-a3) 2cos(k-a2) 

(4) 

where ai==(-l/2, -%/3/2), a2 = (l,0), and 
a3=(— 1/2, \/3/2) represent the displacements in a 
unit cell from A to B site, from B to C site, and from 
C to A site, respectively. In this notation, the first 
BZ is a hexagon with the corners of K=± (27r/3) ai, 
±(27r/3)a2, ±(27r/3)a3. 

The energy spectrum for spinless Hamiltonian i?o(k) is 

characterized by one dispersionless flat band (e[^''=— 2), 
which reflects the fact that the 2D kagome lattice is a line 
graph of the honeycomb structure'*'^, and two dispersive 
bands, e^°^^^ = ItV^ 



3 with fok=X]f=i cos^ (k • ai). 
These two dispersive bands touch at the corners (K- 
points) of the BZ and exhibit Dirac-type energy spectra, 



M 

^2(3)k" 



= (l=F-\/3|k— K|), which implies a "particle-hole" 



symmetry with respect to the Fermi energy €f=1. 
corresponding eigenstates of i?o(k) are given by 



^ik) = C!nk (gik, g2k, gsk)^ : 



The 



(5) 



where the expressions of the components g^k and the nor- 
malized factor G'„(k) for each band are given in Table I. 
At two equivalent BZ edge points M=(0, ±7r/\/3), one 

can find that the wave function |u|jk) is iU defined since 
both its denominator and numerator are zero at these 
two points. 

When an external Rashba SO coupling, which can be 
realized by a perpendicular electric field or by interaction 
with a substrate, is taken into account in the 2D kagome 
lattice model, the spin degeneracy will be lifted. The 
tight-binding expression for this external Rashba term 
can be given as follows 



.A 



{ij)af) 



(6) 



where A is the Rashba coefficient, a are the Pauli ma- 
trices and dy is a vector along the bond the electron 
traverses going from site j to i. Taking the Fourier trans- 
form [Eq. ([3])] and considering the V'k below Eq. ([2]), we 
have i?so=I]k V'k-ffso(k)'(/'k with 



i^so(k) 



(7) 



and 



e'?" sin(k-ai) — e^*"?" sin(k-a3) 

i/_R(k) = A I e*^ sin(k-ai) — «sin(k-a2) 

-e~*t sin(k-a3) — isin(k-a2) 

(8) 

Inclusion of the Rashba SO term in the Hamiltonian 
makes the analytical derivation of the eigenstates \unk) 
(n=l,...,6) and eigenenergies e„k very tedious. At the 
general k points, these quantities can only be numeri- 
cally obtained. At some high-symmetry k points, how- 
ever, they can be approximately obtained, which turns 
out to provide a great help in analyzing SHC. 

The energy spectrum for the total Hamiltonian 
i?(k)=iJo(k)-|-i?so(k) is numerically calculated and 
shown in Fig. 2 (solid curves) along the high-symmetry 
lines (r^K, K->M, and M^T) in the BZ. The Rashba 
coefficient is chosen to be A=0.1. Note that in this paper, 
we only concern the physically reasonable limit of A<Cio 
(to is chosen to be unity). For comparison we also plot 
in Fig. 2 (dashed curves) the energy spectrum in the ab- 
sence of the Rashba SO coupling (A=0). For the middle 
and upper bands, one can see that the spin degeneracies 
are generally lifted in the BZ with the exception at T 
and M points, at which the energy is still spin degener- 
ate due to time-reversal symmetry. The most prominent 
splitting occurs at the corners (K-points) of the BZ. How- 
ever, this splitting does not change the Dirac-type nature 



of the dispersions around these corners. Also, there still 
exists the contacts at these corners between one middle 
band and one upper band, as seen from Fig. 2. For the 
lowest flat band, on the other hand, it reveals in Fig. 2 
that the Rashba splitting is negligibly small, and there is 
no observable SO effect on this flat band. The two-band 
approximation given below will also indicate this fact. 

TABLE I: The expressions for the coefficients in Eq. (5) with 

~^ ^ K^ir-4cos2x2] 

q2k e^°fc cos x-i +2 cos X2 cos X3 

qsk cos X3 + 2 cos X2 cos xi 

G-l 26,eif + [4b, - 3eif ] cos^ x^ + 6(bfe - l)^^ 

The conserved SHC a'^y includes two components, 
^xy—^xy'^'^xy^ whcrc cr*" is the conventional part and 
(T*J comes from the spin torque dipole correction. In 
terms of the band energies e„k and states |w„k) of 

(k)=iJo(k)-l-i?so(k), these two SHC components are 
given bji^ 

< = -eri [/(e„k) - /(e„'k)] 

Im(wnk|^ {vx, Sz} I 

(enk - En'k)^ + rf 

and 

f^Z) = -efi lim — V [/(e„k) - /(e„'k+q)] 

q^o <ix r-f . 

n^n' ,k 

Re(Hnk|T (k, q) |u„'k+q) (Mn'k+q|% (k, q) |u„k) 

(e«k - en'k+q) + f? 

(10) 

where r (k, q) = ^ [r (k) +t (k + q)] with r (k) =8^ , 
V (k, q) is given in the same manner, and f{enk) is the 
equilibrium Fermi function. The limit of 77 — > should 
be taken at the last step of the calculation. In the 
present six-band model the spin operator Sz should be 
written as Isxs^f^ in unit of h/2. 

We have numerically calculated the SHC as a function 
of the electron Fermi energy ep. The main results for 
zero temperature are shown in Fig. 3, in which Fig. 3(a) 
plots the conserved SHC a^y, while Fig. 3(b) plots its two 
components, i.e., the conventional term cr^^ and the spin 
torque dipole term a^y. For comparison, the value of the 
Rashba SO coefficient A used in Fig. 3 is the same as in 
Fig. 2 (solid curves). One can see that within the whole 
range of the electron filling (Fermi energy), the two com- 
ponents a^y and a^y always oppose each other. In fact, 
this feature of opposite signs of the two components crj° 
and a^y (if both of them are nonzero) is robust and does 
not depend on specific models'*^. Remarkably, the ampli- 
tude of a^y is as twice large as that of crj" which results 
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FIG. 2: (Color online) (a) Energy spectrum of the 2D kagome 
lattice with Rashba SO constant A=0.1 (solid curves), (b) and 
(c) show the Fermi surfaces in the regimes — 2<ei?<0 and 
0<eF<l, respectively. Directions of the electron's velocity 
and spin polarization are also shown by red and green arrows 
respectively, (d) Reconstructed Fermi surface around the two 
K points by gluing the six sheets of the Fermi surface in (c). 
For comparison, the energy spectrum in the absence of the 
Rashba SO coupling is also plotted, see the dashed curves in 
(a). One can see that the lowest flat band is immune to the 
Rashba SO coupling. 



in the consequence that the total SHC a^y has an overall 
sign change with respect to the conventional SHC (t*°. 
Together with the previous studies of the conserved SHC 
in the Rashba 2DEG^. As will be shown below, around 
the r point the present model can be mapped into the 
simple Rashba 2DEG model. Thus, one can see the key 
role played by the spin-torque-dipole term, which in some 
special cases tends to overwhelm the conventional SHC 
by an opposite contribution. On the other hand, consid- 
ering the variation of the SHC as a function of electron 
Fermi energy, the present results in the 2D kagome lattice 
display more profound features compared to those in the 
2DEG system. In fact, it reveals in Fig. 3 that the con- 
served SHC and its two components display four plateaus 
as a function ep- When the electron filling satisfies the 
condition — 2.0<ei?<0, the value of a^y is e/Sn, while the 
values of (7^° and cr^y are —e/8n and e/An, respectively. 
When the electron filling increases to satisfy 0<eF<1.0, 
then the conserved SHC jumps down to o'^j,= — e/47r, 
while its two components also jump to cr^^=e/47r and 
a^y=—e/2TT. When the Fermi energy continues to in- 
crease to satisfy 1.0<eF<2.0, then the conserved SHC 
jumps up to a^y=e/4:TT, while its two components also 
jump to a^y— — e/47r and cr*^=e/27r. Finally, when the 
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FIG. 3; (Color online) (a) The conserved SHC a°y and (b) 
its two components a'^y (red curve) and a^y (blue curve) as 
a function of the electron Fermi energy for the Rashba coeffi- 
cient A=0.1. 



Fermi energy satisfies the condition 2.0<ei?<4.0, then 
the conserved SHC jumps down to cr^y= —e/Sn, while 
its two components jump to cr*°=e/87r and cr^^=— e/47r. 

We turn now to understand the physics embodied in 
Fig. 3. Since we are dealing with the usual case of 



weak SO coupling (A^io), thus the SHC behavior in 
Fig. 3 should be mainly due to the coupling of the two 
Rashba SO-split bands and can be described by an ef- 
fective two-band approximation. To be more clear, let 
us treat the Rashba SO term as a perturbation to the 
spinless Hamiltonian _ffo(k). The expressions for the un- 
perturbed eigenenergies e^^^ (n=l,2,3) and eigenstates 
l'"ik) liave been given above. Then, the effective two- 
band Hamiltonian originating from e^*]] and Im^*]]) is ob- 
tained by taking into account the Rashba SO splitting as 
follows 



F„(k) = eS^l2x2 + 



A„ke 







(11) 



where the basis set to expand ffn(k) consists of lu^^J)®] t 
J,). Here the coefficients A„k and ipnk are 



and 



,(o)\ 



defined by 



J 



A„k cos (pnk = 



V3A 



(k)(6i°^+2)(e|^°^'-4cos2 fc,) cos fc, sin(V3fc^), 



(12) 



A„k sin (^„k=-^G2 (k)(e^°J+2) sin [4e^"J cos k,^ + ie^°^^+4. cos^ fc^) cos{V3ky) 

r 



The eigenenergies of Hn (k) are 

(±) _ (0) ^ ^ 



The corresponding eigenstates are given by 



As a result, the total Hamiltonian can now be approxi- 
mated by 



F(k) = ®Li^n(k). 



bands, depending on the position of the Fermi energy. 
Remarkably, there is a particle-hole symmetry between 

(13) the middle and upper bands with respect to their con- 
tact energy plane. As a consequence, the SHC is anti- 
symmetric with respect to the Fermi energy ei?=1.0, as 
revealed in Fig. 3. Keeping this fact in mind, our re- 

(14) maining discussion of Fig. 3 will focus on the two SHC 
plateaus and the transition between them when scanning 
Ei? through the middle bands. According to Eqs. (ISjl- fTU]) 
and our two-band approximation (|lip . when the Fermi 

energy crosses the two middle bands e2k^ it can be shown 

(15) that the conventional part a^^ and the spin-torque-dipole 



This two-band approximation proves to work very well 
in the weak Rashba SO coupling limit. In particular, one 
can see that the lowest flat band (e^k^^ —2) is not split by 
the Rashba SO coupling in the first order in A, since the 
quantity Aike"''!'' is zero and as a result, the off-diagonal 
element in Eq. (|lip {n=l) vanishes. This perturbative 
analysis agrees well with the exact numerical result in 
Fig. 2, which shows that the original flat band e[^'' keeps 
nearly dispersionless upon weak Rashba SO interaction. 
As a result, the contribution of these two spin almost- 
degenerate flat bands to the SHC should be negligibly 
small, which has been verified by our numerical test. 

Thus, the finite SHC in Fig. 3 is ascribed to the contri- 
butions from the two (SO-split) middle or the two upper 



part a^'L of the conserved SHC are given by 



(+) f,Ao) 



2k 



A. 



2k 



(16) 



and 



_^ \ /2k ■/2k 



(0) 
2k 




(17) 



where /2k'' are the Fermi distribution functions for the 
middle bands e2k''- 
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According to the Kubo formulae (fT6l) - ((T7|) . now let us 
see the first SHC plateau in Fig. 3 for — 2.0<ei?<0. Since 
this plateau occurs upon occupation of the bottom (at 
the r point) of the middle bands, thus we can simplify 
the discussion of the first SHC plateau by expanding the 
middle-band Hamiltonian -ff2(k) around the T point up 
to the first order in the Rashba coefficient A 



-2.0 + + X {ky(7x 



kx O'y ) 



(18) 



Not surprisingly, the effective Rashba Hamiltonian ([T| 
around the T point in the present kagome lattice is sim- 
ilar to that in the semiconductor 2DEG. Thus, as has 
been done in the 2DEG system'*^, a straightforward an- 
alytical calculation in terms of Eqs. p6)l - (ll8p gives the 
zero-temperature SHC as cr*"= —e/Sir, o-^^=e/47r, and 
subsequently a^y=e/8Tr. This analytical result is con- 
sistent with the numerical result in Fig. 3 for the first 
SHC plateau. Actually, the first SHC plateau in Fig. 
3 goes beyond this analytical treatment around the F 
point and persists with increasing the Fermi energy up to 
eF=0. The reason is attributed to the equivalent Fermi- 
surface topologies when changing ep within the interval 
[—2.0,0]. In fact, when ep he in the region [—2.0,0], the 
corresponding 2D Fermi surface consists of of two simple 
closed loops circling around the T point, as illustrated in 
Fig. 2(b). Here, from Fig. 2(a) one can see that the crit- 
ical value eF=Q corresponds to the case that the Fermi 
surface nested in the middle bands touches the BZ edge 
at the M point, at which the energies of the two middle 
band are degenerate due to the time-reversal symmetry. 

When the Fermi level goes over this critical value, i.e., 
eF>0, then the Fermi surface abruptly changes its topol- 
ogy. Instead of simple closed loops, the Fermi surface for 
0<eF<l-0 is characterized by six pieces of disconnected 
segments around six corners (K points) of the BZ as 
shown in Fig. 2(c). After gluing these segments together 
by a simple translation operation in the extended BZ, 
which does not change the property of electron states, 
then one can get two sets of closed loops around two K 
points as shown in Fig. 2(d). Thus the number of Fermi 
loops is doubled in the case of 0<eF<l-0 compared to 
the case of — 2.0<ei?<0. This fundamental change in the 
Fermi-surface topology by increasing the electron filling, 
together with the combined fact that (i) the contributions 
from these two sets of K-centered Fermi loops are equiv- 
alent, and (ii) the normal direction of the Fermi surface 
for 0<e_F<1.0 is opposite to that for — 2.0<ei?<0, result 
in a downward jump of SHC plateau from a^y—e/Sir to 
(Jxy= —e/ATT at the critical value of ep—Q. To be more 
clear and to verify this argument based on the Fermi- 
surface topology, near each corner of the BZ let us ex- 
pand the middle-band Hamiltonian H2 (k) up to the first 
order in the Rashba coefficient A, 

/3 

H^ = l-V3k- A^ [kyax - kxUy) , (19) 

where the wave vector k is coordinated with respect 
to the K point. By substitution of the eigenenergies 



and eigenstates of into Eqs. ((T6|) - ((T7)l . and taking 
into account the six corners of the BZ, it is straightfor- 
ward to obtain the zero-temperature SHC as cr*°=e/47r, 
fTxy— ~c/27'', and cTxy— — e/47r, which is consistent with 
the numerical result in Fig. 3. 

Therefore, it becomes now clear that the different 
SHC plateaus in Fig. 3 are due to the different Fermi- 
surface topologies when varying ep- This observation 
makes it highly interesting to reinterpret the metallic 
SHE, like what has been done in discussing the metal- 
lic AHB^i^i^i^, in terms of Berry phases accumulated 
by adiabatic motion of electrons on the Fermi surface. 
The previous work have shown the relationship between 
the SHC and the Berry phase in the Rashba 2DEG^iia. 
The Fermi surface involved in those discussions is as sim- 
ple as shown in Fig. 2(b). Compared to the Rashba 
2DEG, one can see from the above discussions that the 
present kagome lattice provides more profound Fermi- 
surface topologies in the different regions of the electron 
filling. On one hand, in the regime — 2.0<ei;'<0 the effec- 
tive 'T-valley" Hamiltonian ^TE\\ and the Fermi surface 
of the kagome lattice are identical to those of the Rashba 
2DEG. As a result, the two kinds of systems have the 
same Berry-phase SHC in this regime. On the other 
hand, in the regime 0<ei?<l-0 the effective "K-valley" 
Hamiltonian of the kagome lattice, which is absent 
in the Rashba 2DEG, has a remarkable Dirac-type spec- 
trum with linear dependence of the energy on the elec- 
tron momentum. Exploring the K-valley spintronics as- 
sociated with Berry phases is the task of our following 
discussions. 

The Berry phases of Bloch states [u^^^'') for closed paths 
Ck^'^ in the fc-space are written as 



In 



(20) 



where are the Fermi loops identified by the zero- 

temperature Fermi distribution function @{ep — e^"^^), 
and 



-^nk ^ 



'■nk 



(-*9k) 



(21) 



are the Berry connections. The corresponding Berry cur- 
vatures are defined as k'x -^[^ ■ By substituting 
Eq. dUl) into Eq. ©, and noting that afy = -crf^, we 
have 



'xy 



eh 

y 



p.=+- k 
AO) 



Jnk 



nk 



''nk 



''nk ^-^nk 



(22) 



where '^nk~'k~^it' band velocity in the absence of 

the Rashba SO coupling. Now we focus our attention 
to the regime 0<eir<1.0, within which the gluing Fermi 
surface consists of two set of loops around two K points 
as shown in Fig. 2(d). According to the K-valley Hamil- 
tonian (|19p and its eigenenergies e'^ 
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and eigenstates ^2^^^ = ^ (=F*e ^v>2k^i'^^ with (p2k= 

taii~^{ky/kx), it is straightforward to obtain the zero- 
temperature conventional SHC as 



sO_o 



E 



(23) 



where 82'*'' (/i=+,— ) in Eq. ([^5]) denotes the integral 

area bounded by the Fermi loops Cg'*'' [see Fig. 2(d)], 
and the Berry connections 



^2k 



2 dk 



k^' 



k-2k 



(24) 



are equivalent for the two middle bands. Note that the 
factor 2 in Eq. ([221) is due to the contributions from the 
two K valleys. Clearly, if we define an Abelian spin gauge 
field 62k=(0,0,B) withB=[|xA2k]^, then Eq. ^ de- 
notes a spin-flux difference through two areas 82^'' and 
82 ^. In virtue of this way, we define a spin gauge po- 
tential ^2k to satisfy VkX^2k='B2k, then the expression 
for the conventional SHC is rewritten as 



47r2A 



V /i / 62k-dS (25) 
V M / VkxAk-dS 

E^^ <p ^2k-dk. 

We choose a symmetric form for the spin gauge potential 

^2k, 



47r2A — js 



47r2A 



(26) 



which obviously satisfies VkX^2k='82k- By substitution 
of Eq. (Uni) in Eq. dH]), we have 



where /c^^ are the Fermi wave vectors for the two middle 
bands = l—\/Zk ± '\/3A/2, and we have used the fact 



kp 



■kp=\. Thus, we get a remarkable relationship be- 



tween the conventional SHC and Berry phases for the 
K-valley Hamiltonian. Using the chosen middle-band 
eigenstates given above Eq. (|23p . it is simple to obtain 

the Berry phases as 72^''=72 ''^t'', leading Eq. ([77]) to 
o'^?,=i^, consistent again with the numerical result in 
Fig. 3(b). 



In summary, we have theoretically investigated the 
metallic spin-Hall effect in the 2D kagome lattice with 
Rashba SO coupling. When varying the Fermi energy 
ef, we have found that the conserved SHC cr^y and its 
two components, i.e., the conventional term cr^^ and the 
spin-torque-dipole term cr*^, are characterized by a series 
of plateaus, which is absent in the simple 2DEG system. 
In the whole range ep varies, the two terms cr'^y and (7^]^ 
have opposite contributions. The magnitude of cr^]^ is 
twice of that of cr^° . It has been shown that these SHC 
plateaus in the different regions of are closely associ- 
ated with the topologically different Fermi surfaces sur- 
rounding the high-symmetry BZ points, i.e., the T and 
K points. Thus, as has been revealed in this paper, a re- 
lationship between these SHC plateaus and Berry phases 
accumulated by adiabatic motion of quasiparticles on the 
Fermi surfaces can be built up, which is similar to the 
metallic AHE. In particular, we have shown that com- 
pared to the case with the Fermi surface encircling the T 
point, the amplitude of the SHC with the Fermi surface 
encircling the K points is twice as large. Considering the 
combined fact that (i) the 2D kagome lattice is line graph 
of the honeycomb structure, (ii) the Rahsba SO coupling 
and the Fermi surface surrounding the K points can be 
easily realized in the graphene with honeycomb structure, 
and (iii) the similar Berry-phase AHE has been recently 
observed, we expect that the present prediction of the 
K-valley enhanced SHE can be observed in the graphene 
system. 
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